TARGET FRAGMENTATION AND FRACTURE FUNCTIONS. 
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We analyse recent data on the production of forward neutrons in deep inelastic 
scattering at HERA in the framework of a perturbative QCD description for semi- 
inclusive processes, which includes fracture functions. 



In the most naive quark-parton model picture, the semi-inclusive cross 
section for the production of a hadron h from the deep inelastic scattering of 
charged leptons is usually taken to be given byld: 
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where, fi/ p is the parton distribution of flavour i, D^u is the fragmentation 
function of a hadron h from a parton i, and Cj = Aire 2 q .a 2 / x(P + I) 2 . The 
kinematical variables used to characterize these processes are: 
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where q is the transferred momentum (— q 2 = Q 2 ), I and P are the incoming 
lepton and proton momenta respectively. Eh, E p and Oh are the produced 
hadron and target nucleon energies, and the angle between the hadron and the 
target in the centre of mass of the virtual photon-proton system, respectively. 

Although next to leading order corrections to this cross section are also well 
known, and have been shown to give a very good description of the so-called 
current fragmentation region {Oh, > t/2), the target fragmentation region, 
which corresponds to Oh — (z^ — 0), cannot be described with this simplified 
picture. First of all, it is easy to see that, at the lowest order, hadrons can 
only be produced antiparallel to the target nucleon (Oh — t), excluding the 
forward configurations. On the other hand, going to next to leading order, 
the corrections to the cross section develop divergences proportional to l/zh, 
related to soft emission (Eh = 0), and also to collinear configurations where 
hadrons are produced in the direction of the remnant target (Oh = 0). Since at 
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lowest order hadrons cannot be produced in that direction, it is not possible to 
factorize the divergence, as usual, into parton distributions and fragmentation 
functions. Then, in order to describe hadrons produced in the target frag- 
mentation region even at the lowest order, and also to be able to perform at 
higher orders a consistent factorization of divergences originated in the current 
fragmentation region, a new distribution_.has to be introduced, the so-called 
fracture functions, M ith / N (x,(l — x)z)aR These distributions represent the 
probability of finding a parton of flavour i and a hadron h in the target TV 
(here z = E h /E p (l - x)). 

Therefore the complete leading order expression for the cross section be- 
comes 
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Higher order corrections to this kind of cross section can be found in refs. ifii. 

The scale dependence of fracture functions at 0(a s ) is driven by two kinds 
of processes, which contribute to the production of hadrons in the remnant 
target direction: the emission of collinear partons from those found in the 
target (the usual source of scale dependence of parton distributions, often 
called homogeneous evolution), and those where partons radiated from the 
one to be struck by the virtual probe, fragment into the measured hadron (the 
so-called inhomogeneous term). These two contributions lead at leading order 
to the following equation: 
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where P]{u) and Pj' l (u) are the regularized and real Altarelli-Parisi splitting 
functions, respectively. 

Recently, the ZEUS CollaboratiorLhas measured DIS events identifying 
high-energy neutrons in the final statelj, at very small angles with respect to 
the proton direction (9i a b < 0.75 mrad), in the kinematical range given by 
3 x 1(T 4 <x<6x 1CT 3 , 10 < Q 2 < 100 GeV 2 and high x L ~ z(l - x) > 0.30. 
The ZEUS Collaboration have reported that events with xl > 0.50 represent 
a substantial fraction (of the order of 10%) of DIS events. In the framework of 
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a picture for semi-inclusive processes including fracture functions, as the one 
outlined above, the ZEUS findings can be represented (at LO) by 
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In eq.(5) we have also defined the equivalent to F2 for fracture functions: 
M^ p (x, xl, Q 2 ) = xJ2i e2 Mi.n/p( x , x L,Q 2 ) and we have made explicit the 
integration over a finite (measured) range of x^. 

In fig. 1 we show the experimental outcome for this fracture function (as 
defined in eq. (5) and at Q 2 — 10 GeV 2 , taking advantage of the negligible 
Q 2 dependence of the data), and we compare it to i 7 ^ and F£. We also show 
the contribution to the same observable coming from current fragmentation 
processes, which is pure NLO and it is about 8 orders of magnitude smaller 
than the experimental data. 

Fracture functions, as parton distributions in general, are essentially of a 
non-perturbative nature and have to be extracted from experiment. However, 
their close relation with fragmentation and structure functions allows in certain 
extreme cases a model estimate for them, which can then be compared with 
actual measurements and evolved with the corresponding evolution equations. 
As an example, recently,_a very sensible model estimate for the production of 
forward hadrons in DIS 0, exploiting the idea of non-perturbative Fock com- 
ponents of the nucleon has been proposed. In this approach the semi-inclusive 
DIS cross sections, and through them the corresponding fracture functions at a 
certain input scale Qq, can be interpreted as the product of a flux of neutrons 
in the proton (integrated over p^) times the structure function of the pion 
exchanged between them, i.e. 

M 2 re/P (x, x L ,Q 2 ) ~ n/p (x L ) Ff (y^, Qo) • (6) 

Using a non-perturbative computation of the flux, which is in very good agree- 
ment with experimental data on high energy neutron and A ++ production in 
hadron-hadron collisions and a parametrization for the pion structure function 
□, in fig. 1 we make a comparison between the model estimate and the data 
for Qq = 10 GeV 2 , finding a remarkable agreement between them. 

The success of the model estimate encourages us to go further and use 
the functional dependence of fracture functions, induced by the model and 
corroborated by the data, to analyse also the Q 2 dependence. 

We first analyse the more familar process of neutron production. Since 
the probability of current parton fragmentation into a neutron (given by frag- 
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Figure 1: Fracture function of neutrons in protons as measured by ZEUS compared to the 
model prediction and current fragmentation contributions 



mentation functions) is comparatively small with respect to that of processes 
originated in the target (fracture functions) , no significant effects are expected 
in the scale evolution arising from the inhomogeneous term in this process. 
The evolution is mainly driven by the usual homogeneous term of the evo- 
lution equations leading to an almost constant ratio between the number of 
neutron tagged events and that of all DIS events, as observed by ZEUS. 

However, the scale dependence induced in the cross section for the pro- 
duction of pions, at least in the kinematical region of very small x and small 
xl, can be considerably affected by the inhomogeneity, given that soft pions 
are produced more copiously from quarks than from neutrons. In order to 
analyse these features of the evolution, wc estimate the proton to pion frac- 
ture function at some input scale Qq using the same ideas formerly applied to 
neutron production, and noticing that the flux can be straightforwardly ob- 
tained from the one used in the last section by means of the crossing relation 
4>tt+ /p( x l) = 4>n/p(^ ~ x l)- Then, the proton to pion fracture function can be 
approximated by 
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In fig. 2a we show the model estimates for proton to pion fracture functions 
(taking Qq — 4 GeV 2 ), integrated over two different bins of xl, compared 
with the contribution coming from the current fragmentation processes. We 



assume here the same restrictions as in the data from the ZEUS Collaboration 
for neutron production. 




Figure 2: a) Prediction for the fracture function of n + in protons for two different bins of 
ij, and the current fragmentation contribution, b) Evolution of the fracture function of 7r + 
in protons for x L = 0.50 and Q 2 , = 1 GeV 2 , c) x L = 0.10, and d) Q\ = 4 GeV 2 

Of course, the model is not expected to work over the whole kinematical 
region and, in fact, any deviation from the scale dependence implied in eq. 
(8) (note that the flux is assumed to be Q 2 independent) would show the 
breakdown of the approximated factorization hypotesis. However, the ansatz 
in eq. (8) can be taken as an effective relation, valid at some initial value of 
Qq, for which the estimated flux is adequate, and therefore provides a sensible 
input distribution. As usual, the correct scale dependence is that given by the 
evolution equations for fracture functions, and that is the aim of our next step. 

In order to study the effect of the inhomogeneity in the evolution we take 
different values of Xl, and keep them fixed while we analyse the x and Q 2 
dependence of fracture functions induced by both the homogeneous and the 
complete evolution equations. In fig. 2b we show the result of an evolution 
from Q 2 — 1 GeV 2 a,t xl — 0.50. Both solutions, the homogeneous (dotted 
line) and the complete (solid line) are superimposed, the difference being less 
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than 0.1%. These behaviours are perfectly compatible with the results ob- 
tained by the ZEUS Collaboration in the case of neutron production (where 
the inhomogeneity contributes about 10 times less) in the same kinematical 
region, where no difference has been found in the evolution between F2 and 
M 2 . 

However, for smaller xl the situation is completely different. As the frag- 
mentation function increases with lower values of the argument, the inhomoge- 
neous contribution becomes much more relevant and its effect in the evolution 
is sizeable. In fact, fig. 2c shows the evolution result for xl — 0.10 and Qq = 1 
GeV 2 , where the full evolution results outsize the homogeneous one by a factor 
of 4 at small x. These corrections are smaller if the ansatz of eq. (8) is assumed 
to be valid at values of Qq = 4 GeV 2 (fig. 2d) but still remain considerable. 

Concluding, in this paper we have analysed recent experimental data on the 
production of forward neutrons in DIS in terms of fracture functions, finding 
that the main features of the data can be fairly reproduced by this perturbative 
QCD approach, once a non-perturbative model estimate for the input fracture 
functions is given. Studying the evolution properties of these fracture functions 
in the specific case of forward pions in the final state, we have found that the 
effects of the inhomogeneous term in the evolution equations are large and 
measurable, particularly in the kinematical region of very small x and small 
xl. These effects are negligible for large values of Xl, justifying the use of 
the usual homogeneous Altarelli-Parisi equations for, as an example, the t- 
integrated diffractive structure function F®^ {x^ , /?, Q 2 ), which is just the 
fracture function of protons in protons M^ v (fix^, (1 — x ip)i Q 2 )i. 
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